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F/A-18 HARV PID REPORT 


INTRODUCTION 

The following is a compilation of PID results for both longitudinal and lateral directional 
analysis that was completed at West Virginia University during the Fall 1994 semester. 

During Summer 94 it was established that the maneuvers available for PID containing 
independent control surface inputs from OBES were not well suited for extracting the 
cross-coupling static (ie. Cj^j ) or dynamic (ie. Cj^pf) derivatives. This was due to the fact that 
these maneuvers were designed with the goal of minimizing any lateral directional motion during 
longitudinal maneuvers and vice-versa. This allows for greater simplification in the aerodynamic 
model as far as coupling between longitudinal and lateral directions is concerned. As a result, 
efforts were made to reanalyze this data and extract static and dynamic derivatives for the 
F/A-18 HARV without the inclusion of the cross-coupling terms such that more accurate 
estimates of classical model terms could be acquired. 


FLIGHT DATA ANALYZED 

The following listings show the flight numbers, maneuver types and angle of attack 
ranges examined in this analysis. 


Longitudinal: 


Flight 

Maneuver Type Angle of Attack (deg) 

166 

SSI MD 

25,30,40,50,60 

226 

SSI MD 

20,30,40,50 

250 

OLON 5 e =7.5° 

20,40 


3211 $ e =7.5° 

20 


1" STICK PITCH STEPS 

20,30 

253 

SSI MD 

20,30,40,50,60 

Lateral: 

Flight 

Maneuver Type Angle of Attack (deg) 

226 

SSI MD 

10,20,30 

248 

SSI MD 

10,20,30,50 

250 

SSI MD 

50,60 


OLAT 5 a =5°,5 r =10° 

20,40 


♦SSI MD - Single Surface Input Multiple Doublet maneuvers (OBES ON, RFCS ON) 

*OLON 5p=7.5° - Optimal Longitudinal Input, elevator deflected 7.5° (OBES ON, RFCS ON) 
*3211 5 e = 7.5° - 3211 Frequency Sweep, elevator deflected 7.5° (OBES ON, RFCS ON) 
*1" STICK PITCH STEPS - Pilot applied 1" control stick pitch steps (OBES OFF, RFCS OFF) 
*OLAT 5 a =5°,5 r = 10° - Optimal Lateral Input, aileron deflected 5°, rudder deflected 10° 

(OBES ON, RFCS ON) 



Note that each flight often has more than one run at each angle of attack tested. The 
proper time cuts, maneuver descriptions and average values of all required mass characteristics 
(lx, Iy, etc. ) for each time cut have been placed in a textfile database available upon request. 


LONGITUDINAL SETUP 

Four longitudinal flights containing static PID maneuvers were examined. Time histories 
of leading-edge and trailing-edge flap deflections clearly indicated high correlations between 
themselves, as well as angle of attack, due to their scheduled deflections based on a. In order 
to avoid encountering any unwanted biases in the estimated derivatives all leading-edge flap 
derivatives (C N51ef , C m51ef & C A51ef ) were held constant at the appropriate wind tunnel values. 
All trailing-edge flap derivatives (C N5tef , C m5tef & C Agtef ) were also held constant at 
appropriate wind tunnel data ONLY if the trailing-edge flaps were NOT PULSED by OBES 
If die trailing-edge flaps were independently pulsed by OBES then these coefficients were also 
estimated. All wind tunnel values were based upon angle of attack and Mach for the maneuver. 
Also note that all wind tunnel data was presented in terms of CV , Cr) and C . The lift and 
drag terms were easily converted to normal and axial force coefficients by using: 

C N =C L c os (a) +C D sin (a) 


c a =c d cos (#) -Cj.sin(o) 


All mass properties (x 


y C g> z cg> ^x’ ^y ’ ^7> i 


- - ... eg* >cg’ z cg» A x’ V *Z’ ixz ^ mass ) were considered to be 
unique values for each individual maneuver. These values were also cut using getFdas and an 
average value of each property over the maneuver time was used in pEst for the model. 

The only dynamic derivatives estimated here were C Nqf , C m f and C Aaf due to the lack 
of excitation in the lateral direction during the maneuvers. The classical Static derivatives 
estimatedwere. C No , C^ a , C^, c N5tef> C N5pv^ C NSyv’ c N5as> c mo’ C ma- c mSe’ 

mStef’ m5pv’ ^m5yv’ ’-mSas’ c Ap’ c Aa> c A5e’ C A6tef’ C A6pv’ C A5yv’ C A5as 

Two additional parameters anBias and axBias were estimated to model the measurement 

biases for the normal and axial linear accelerometers. The goal with respect to the longitudinal 
cost function (J) was to minimize the difference between measured and computed values of: a, 
<1> a* and climb angle y. Note that although q mo< j was computed no weighting was 

applied in the cost function for it. This was done because q and q mfv) are nearly identical 
responses where: 1UU 


< 3 m od =< 7- w sin(p) 

In these longitudinal directional OBES maneuvers o> and j8 are near zero. Consequently, if q 
measured and q computed match properly then naturally q mod will follow by definition. 


LONGITUDINAL MODEL 

The classical state equations already available in pEst for alphadot, qdot and thetadot 
were used. Note that the new state equation gammadot describing directional change in the 



velocity vector (climb angle) was added. It is a function of the total coefficients Cj^, C A and 
Cy. 

Y = -2L [n x (sinacosii+coaasinPsini») -n y (cospsinti) +n x (cosacosii-sinasinPsin|i) -cosy] 
V T 


where: 


x=- — C A + — n y= — C Y n z= — 

x mg A mg y mg mg 


While decoupling longitudinal and lateral, as was the goal with these maneuvers, it was found 
that the Cy term could be safely set to zero and thus its subcomponents not estimated. This 
approximation holds since n, the bank angle about the velocity vector, is often small during such 
maneuvers. The total coefficient buildup for the longitudinal model can be found in [2,3] while 
the added state equation for 7 developed by Kalviste using spherical mapping, may be seen in 
[1,3]. 

Also recall that the model being used breaks aircraft rotation down into three of the four 
components possible at each discrete time point (p mcK j* q mo d> r mod allowing for the 

estimation of coefficients in the buildup equivalent to rotary balance (ie. Cj w , etc.) and forced 
oscillation derivatives (ie.C N qf, etc.). Details on this modeling technique can be found in [2,3]. 


LATERAL SETUP 

Three lateral directional flights of PID static maneuvers were also examined. In this case 
it was found necessary to hold the differential horizontal stabilizer derivatives constant to wind 
tunnel data if this control surface was NOT pulsed independently by OBES. This was found to 
be the case in the "optimal lateral input" runs made using OBES where some correlation, 
although small, was still present between aileron pulses (made using OBES) and differential 
stabilizer. As before, the wind tunnel data was based on Mach and angle-of-attack for the 
considered run. 

All mass properties were again taken to be unique values for each individual maneuver. 
These values were also cut using getFdas and an average value of each property over the 
maneuver time was used in pEst for the model. 

The dynamic derivatives estimated were Cy^f, Cy^, Cy w , C^f, Cj^, Cj w , C np f, C nr y 
and C nw due to the lack of deliberate longitudinal excitation in the lateral maneuvers. The 
classical lateral static derivatives estimated were: Cy 0 , Cyo, Cyg A , Cy^jj, Cy^, Cyg pv , 

C Y5yv> C kv C 15A> C 15dH> C 15R> C 15pv> C 15yv> C no> ^n/3* C nSA’ C n5dH> C nSR’ C n<$pv’ 

The additional parameter ayBias was estimated to model the measurement bias for the 
lateral directional accelerometer. 

The goal regarding the lateral cost function (J) was to minimize the difference between 
measured and computed values of: /3, p, r, 0 , ay, p m<x j» r mod» 03 and As OPP 0 ^ t0 the 
longitudinal case, the time history for p is not similar to p m(X j and r is not similar to r m(X j. 
Thus, each of these responses had to be included in the cost function. 



LATERAL MODEL 

The classical state equations already available in pEst for betadot, pdot, rdot and phidot 
were used. The new state equation mudot describing bank angle 00 of the aircraft about the 
velocity vector was added. It is a function of the total coefficients C A and Cj^, body axis roll 
rate (P), body axis yaw rate (R) and the rate change of heading for the velocity vector 
(sigmadot). Due to the design of the available lateral directional maneuvers any deviation from 
the longitudinal steady state trim is minimized. Since this proposes a decoupling of longitudinal 
and lateral motion the C A and C N terms are approximated as zero (their subcomponents are not 
estimated). Examination of the relation for mudot shows that the contribution due to these two 
terms are small compared to those of P and R. 

|i=(P coso+P sina) secP+*siny + -^ (n^ino+rijCosa-cosiicosY) tanp 

V T 

In this analysis sigmadot was set in pEst with the option to be computed as a response. 
Sigmadot is a function of the total force coefficients C A , Cy and Cj^. Note that sigmadot was 
not integrated with the goal of matching a measured and a computed because it is equivalent to 
a compass heading. This would be an awkward response to match. Consequently, it made more 
sense to try matching sigmadot measured and sigmadot computed (or the rate change of the 
velocity vector heading). However, here it turns out that C A and Cj^ have a strong contribution 
in the model for sigmadot as seen below. 

ft- 3. [n (sinasinii-co8«sinPcos|i) +n Y (cospcos|») +n,(coB«sinii+sinasinpcosii) ] 

Vycosy x 

As a result, the assumption of a decoupled model where C A and Cjq have negligible contribution 
is not appropriate. Consequently, the measured time history for sigmadot was used in the model 
throughout the analysis. Recall the previous expressions shown for n x , ny and n z and again note 
that the lateral directional aerodynamic buildup expressions can be seen in [2,3] while the two 
new state equations developed by Kalviste can be seen in [1,3]. 


CONCLUSIONS 

A more detailed evaluation of the results are planned to be incorporated into a proposed 
journal paper during Spring 95. Overall, with the data analyzed it was found that less scatter, 
as well as smaller Cramer-Rao bounds, were achieved for the classical parameters over the 
previous attempt at modeling the cross-coupling derivatives in strictly longitudinal and lateral 
maneuvers. In addition, the majority of the control surface coefficients were well estimated 
especially with the use of the SSI MD inputs. One problem encountered with the independent 
longitudinal pitch vane pulse at the end of the appropriate time cuts was that the system 
experienced difficulties returning the computed model back to trim after recovery from the pulse. 
As a result, the longitudinal SSI MD time histories were adjusted to include elevator, trailing- 
edge flap and symmetric aileron only. Perhaps a pitch vane doublet would prove better than a 
single pulse since the yaw vane doublets did not encounter this problem. 

Enclosed with this document are the full set of longitudinal and lateral-directional 
parameter estimate plots showing coefficient estimates along with Cramer-Rao bounds. All 
bounds were multiplied by a factor of three to account for any modeling errors. In addition, a 
representative time history match for each type of maneuver tested at each angle of attack is also 



enclosed. Additional plots of time histories are available upon request. 
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Abstract 

The subject of this paper is the determination 
of the aerodynamic stability and control derivatives 
from flight data using the Maximum Likelihood 
method for the NASA F/A-18 HARV at high angles 
of attack. The parameter identification (PID) code 
pEst, developed at NASA Dryden, was provided by 
NASA and modified for an alternative modeling 
approach for high angle of attack conditions. 
Estimates were obtained for longitudinal dynamics 
parameters. 

Symbols 

ajj = normal acceleration, g 

= longitudinal acceleration, g 
ay = lateral acceleration, g 

b = wing span, ft 

~c = mean aerodynamic chord 

C- = aerodynamic coefficient where 

i = (L,N,D,A,m,Y,l,n,6), rad -1 or deg -1 
E = expected value 

g = gravity acceleration, ft/sec^ 

J = quadratic cost functional 

= upwash correction factor 
= sidewash correction factor 
= aircraft mass, slugs 
iif = number of discrete time points in a 
maneuver 

n z = number of responses in the cost 
function J 

nt = number of parameters estimated 

P(z) = maximum likelihood probability 

p = roll rate, deg/sec 

q = pitch rate, deg/ sec 

cf = dynamic pressure, lbs/ 

r = yaw rate, deg/sec 


# Assistant Professor, AIAA member 
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wing planform area, ft 
velocity, ft/sec 
variance 

response error weighting matrix 
"a priori" weighting matrix 
X-axis position, ft 

computed aircraft response vector, or 
Y-axis position, ft 

measured aircraft response vector, or 
Z-axis position, ft 

Greek 

angle-of-attack, deg 

sideslip angle , deg 

increment 

deflection, deg 

pitch attitude, deg 

parameter vector to be estimated 

"a priori" parameter estimates 

correlation coefficient 

summation 

roll attitude, deg 

yaw attitude, deg 

angular rate vector, deg/sec 

natural frequency, rad/sec 

gradient 

Subscripts 

steady state trim flight conditions 

axial force 

aileron 

center of gravity 
drag force 
elevator 
stabilizer 
lift force 
rolling moment 
pitching moment 
modified angular rate 
normal force 
yawing moment 



R = rudder 

Y = side force 

Superscripts 
T = transpose 

-1 = inverse 

= first derivative with respect to time 
— = vector quantity 

Introduction 

Today’s high performance military aircraft 
are required to be able to fly at high angles of attack, 
often at high angular velocity. 

The National Aeronautics and Space 
Administration (NASA) is currently involved in a 
high alpha technology program using different aircraft 
which has the following objectives: 

-to provide flight validated prediction methods 
including experimental and computational methods 
that accurately simulate high angle of attack 
aerodynamics, flight dynamics and flying 
qualities 

-to improve aircraft agility at high angles 
of attack 

Part of the overall program for each test 
aircraft involves a parameter identification (PID) 
analysis in which a complete set of stability and 
control parameters is extracted from flight data at 
high angles of attack. A first purpose for this 
analysis is to evaluate the correlation of the flight 
data estimated from aerodynamic characteristics with 
the correspondent wind tunnel and/or numerically 
predicted characteristics. An additional purpose of 
this analysis is to extend the data base of the aircraft 
aerodynamics which, in turn, allows to update the 
simulation codes in flight simulators and provide the 
control engineers with a better understanding of the 
aircraft non-linear behavior for a more efficient 
design of non-linear control laws. 

The use of flight data to estimate aircraft 
stability and control derivatives has been implemented 
for many years. The most efficient parameter 
identification (PID) method to date is the Maximum 
Likelihood (ML) approach used in conjunction with 
a Newton-Raphson (NR) algorithm for the 
minimization of a cost function related to the 
estimation process. This method has been used at 
NASA Dryden since the late 1960’s and several 
codes have been developed for the purpose. 

This paper is relative to the determination of 
the stability and control derivatives for the NASA 
F/A-18 HARV at high angles of attack using the ML 
method with NR algorithm. The NASA F/A-18 


High Alpha Research Vehicle (HARV) is a special 
research aircraft equipped with a thrust vectoring 
system consisting of axisymmetric nozzles and six 
post exit vanes that allow vectoring capability in both 
pitch and yaw. An additional feature of the aircraft 
which makes it very attractive for PID investigation 
is the onboard excitation system (OBES) in the flight 
control computer which allows independent control 
surface motions. The PID process has been 
performed using the pEst software, which is the more 
recent and sophisticated estimation code developed at 
NASA Dryden. Among the capabilities of the pEst 
code there are efficient user interaction and the 
handling of non-linearitites in the dynamic equations 
of motion. 

It is clear that at high angle of attack 
conditions the concept of evaluating the stability 
derivatives about a steady state flight condition can 
no longer be applied. This is due to the inherent 
aerodynamic non-linearities as well as the potential 
coupling between the longitudinal and lateral 
dynamics. The main contribution of this paper is to 
implement a particular aerodynamic and dynamic 
modeling procedure for these non-linear 
conditions^ used in conjunction with the ML 
method for PID purposes. 

This investigation has evolved from a 
previous PID analysis using a more conventional 
modeling procedure. ^ 

Overview of the Maximum Likelihood Method with 
the Newton-Raphson Algorithm 

The most successful parameter identification 
approach to date has been the ML method with a NR 
technique to minimize a quadratic cost functional 
containing differences between the aircraft measured 
and computed responses. This implies the 
maximization of the probability that the aircraft 
model computed responses, using a given set of 
estimates for the unknown aerodynamic coefficients, 
are representative of the true system dynamics. The 
relations forming a basis for this technique are shown 
in Figure 1.^ Define P(z/£) as the probability that 
the actual system response z occurs for a given value 
of the unknown parameters placed in a vector £. 
Define P(£) as the probability that the unknown 
derivatives in the parameter vector match some "a 
priori" values. It is assumed that P(z/£) and P(£) are 
independent and follow Gaussian distributions. 

It is clear that, as the accuracy of the 
estimates increases, the differences between the 
values of the components of z and y, at the same 



discrete time index, decrease. Using the 
independence assumption for the two probabilities 
yields the ML cost functional P(z). The cost function 
J(£) must therefore be minimized in order to 
maximize the ML cost functional. A modified NR 
algorithm is used to solve the associated system of 
equations by using the first and second gradients of 
the cost function with respect to the vector 
The relations governing the NR algorithm are listed 
in Figure 2. The process is iterative with updating of 
the vector £ until some user defined convergence 
criteria is met. Using a Taylor Series expansion an 
expression is generated for J(£). Using the Newton- 
Raphson method, this cost functional is minimized by 
setting its gradient with respect to £ equal to zero. 
This allows for the determination of an iterative 
expression for £ with the first and second gradients of 
J(£). An overview of the aircraft parameter 
estimation process is shown in Figure 3. 

Aerodynamic and Dynamic Modeling at Non-linear 
Conditions 

One of the main purposes of this paper is to 
introduce a particular modeling procedure for the 
non-linear dynamic and aerodynamic conditions of the 
F/A-18 HARV at high angles-of-attack. The 
modeling starts from the classic non-linear 
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longitudinal and lateral-directional equations of 
motion along the stability axes as seen in [ 6 , 12 ]. 
They are shown in Figure 4. A particular innovation 
described in [ 3 , 4 ] consists of the introduction of a 
new set of dynamic axes called the dynamic stability 
axes. A key characteristic of this new dynamic 
frame is that the axes are not orthogonal. By 
definition, the equations of the rotation of the aircraft 
around these dynamic stability axes are given by: 
-Equation of the time rate of change of the flight path 
bank angle (ji ) 

-Equation of the time rate of change of the flight path 
elevation ( 7 ) 

-Equation of the time rate of change of the flight path 
heading angle (a) 

These equations are shown in Figure 5. These 
relations make use of the body axis angular rates to 
obtain the transformation to the dynamic stability axis 
rates. 


Having defined the new system of aircraft 
dynamic equations new aerodynamic coefficients 
describing this motion must be introduced. The 
classical stability parameters (Cjq,C A ,C m ,Cj,C Q and 
C y ) can be separated into two components. Each of 
the total non-dimensional stability derivatives can be 
partitioned into static and dynamic terms: ^ 
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Figure 1: Relations forming a basis for the Maximum Likelihood method 
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Figure 2: Relations governing the parameter update process using the Newton-Raphson algorithm. 
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C l Ja,p,VAt>J>,QJt,a,il) (23) 

where i = (N,A,m,l,n and Y). The static terms are 
made up of non-rotational motion while the dynamic 
terms contain all aircraft motion parameters. The 
static terms are the known conventional derivatives 
often based on aircraft geometry which may be 
represented as follows in their body axis forms: ^ 


Cjt m t 

(24) 

r ^i.+C lt _g*C 4f P»C jl .S 

(25) 

c -~“ c -. +c «. a+c », p+c «. 5 

(26) 

C l =c, +c,p +C,a + c it a 

(27) 


C m -C +C.P+C.B+C.4 (28) 

n 0 Up* Hg K| 


C Y ^C r +C r p +C Y a +Cy* (29) 

Note the cross-coupling stability derivatives 
in these relations. These are required due to the 
presence of dynamic coupling between the 
longitudinal and lateral directions of motion at high 
angles-of-attack. The dynamic terms are new and 
warrant further examination. 

The dynamic derivatives may be built up 
from components due to each rotational term as 
shown in Figure 6.^ Estimates for these terms are 
found during conventional wind tunnel tests utilizing 
the rotary balance and forced oscillation methods in 
addition to translational acceleration terms. The total 
contribution of these dynamic terms are also shown 
in the component build-up of Figure 6. The first 
term in the equation is obtained from rotary balance 
wind tunnel tests. The next three terms are measured 
in the oscillation wind tunnel test. The last two terms 
are the acceleration derivatives. The P m(X p Qxnod 
and R moc j terms are the components of total rotation 
vector about the x, y and z axes. The total rotation 
vector of the aircraft has been divided into only three 
components. The reason for this is that division into 
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(16) 

— 

<j> = P + tan0(Pcos<J> + <2sin<t>) 

(17) 


Where: 



C L =C,f os(o) -C,sin(a) 

(18) 


C c =C A cos(a) +C N sin(a) 

(19) 


Figure 4: The classical non-linear longitudinal and lateral equations of motion about the stability axis. 
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Figure 5: Relations transforming body axis angular rates to dynamic stability axis rates. 
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Figure 6: Relations showing the dynamic axis stability derivative build-up. 


four components along the X, Y, Z and velocity 
vector often results in a modified rotational rate that 
may be greater in magnitude and opposite in sign 
than the measured body axis data. Entering roll rates 
into the equations of motion with the wrong sign 
would cause misinterpretation of the actual aircraft 
dynamics. The method for breaking the total 
rotational vector into only three of the four available 
components result in the following modified roll-rate 
equations: 


=P- wcos acos P 

(32) 

<?«od = <? _USin P 

(33) 

R mod =R- 0 ) sinacosp 

(34) 


The procedure begins by setting the modified 
roll-rate to zero and solving for u>. The 
corresponding values of and Rjj^ are then 

calculated. If they are less in magnitude compared to 
the total body axis values for Q and R, and of the 
same sign, then the rotation vector is broken into 

Qmod> R mo d ^ If test fails then toe same 

procedure is performed with the Q m(X j and R moc | 
relations. Should all three methods fail then the 
rotation vector is broken up into the body axis P, Q 
and R values. The result of these tests at each time 
point are placed into the dynamic parameter model 
previously shown where it affects the overall aircraft 
response. The translational acceleration terms found 
in the dynamic build-up are defined as: 

d r =d -<?+(Fcosa +Rsina)tanp (35) 

(36) 


where the first term in each relation contains effects 
due to aircraft velocity vector translation and body 
axis rotational motion. The trailing terms represent 
effects due to aircraft rotational motion alone. 
Again, the total non-dimensional stability derivatives 
are the sum of the static and dynamic components. 
The total components for Cj, C m and C fl affect the 
aircraft motion through the body axis equations for 
P, Q and R. The terms for C, and C D are located 
in the expression for V while Cjq, and Cy affect 
the body axis response equations for a z , a x and ay 
respectively. 

Modifications to the pEst code to include 
these coupling static and dynamic terms, as well as 
the general equations allowing for aircraft flight path 
angle translation, have been implemented for 
integration with the parameter estimation process. 

Results From the Parameter Estimation Process 

One flight was performed with the NASA 
F/A-18 HARV at six different angle of attack trim 
conditions (10°, 25 \ 30°, 40°, 50\ 60°) for use 
in extracting the longitudinal aircraft stability 
derivatives. The longitudinal derivative components 
of normal force (C N ), aircraft pitching moment (C^) 
and axial force (C^) were chosen to be estimated in 
the PID process. The data relative to the longitudinal 
flights were run through the pEst code to produce 
parameter estimates along with measured and 
computed response curve-fits. The quadratic cost 
functional was set to include differences between 
measured and computed a, q, 0> a^ a x , y and 
q mo d- The response error weighting matrix W ^ was 
an 8x8 diagonal matrix. Response errors for 
measurements with higher measurement noise such as 
a, a and a x were given much lower weightings than 
for those responses recorded with a lower noise level. 
The "a priori" weighting W 2 was not implemented in 
the PID process. 

A very useful feature in pEst is the 
capability to generate the Cramer-Rao bound. This 


P r =p-Psma+Rcosa 



is an experimental variance, presented to the pEst 
user in standard deviation form, which can be used as 
a measure of accuracy for the estimates for the given 
data. For each parameter this bound gives the 
mini mum variance of the estimate using: 


vaiil) 2 



(37) 


where, for a normal distribution: 



(t-P 
„ 2 


(38) 


To account for modeling discrepancies the 
Cramer-Rao bounds are multiplied by a factor of 10. 
Their most practical use is as a measure of the scatter 
of the final estimates for the given set of flight data. 

The results of the PID have shown a 
remarkable agreement between the actual time 
histories and the ones generated with the 
mathematical model using the estimates. A 
comparison of the actual and computed a, 0 and q 

time histories at a = 10°, 25 , 30 , 50 are 

shown in Figures 7-10. The matching between the 
actual and computed time histories for 6 is excellent, 
due to the accurate and noise-free 8 measurements. 
The matching is acceptable for the a and q 
measurements, which show the presence of a 

substantial level of noise. 

The main target derivatives in this 
investigation included C Nq , C^ 5E , and C m5E . 
The results are shown in Figures 11-14. The 
outcome of this investigation has been compared with 
previous results^ obtained with convential modeling 
procedures using the same flight data. Some 

reservations exist on the goodness of this flight data 
at a = 25°. It is believed that at that particular 
condition the aircraft longitudinal dynamics had not 
been properly excited due to some deflection of the 
trailing and leading edge flaps not coordinated with 
the elevator deflections. The results of the PID 
process for Cj^ a shown in Figure 1 1 indicate that the 
introduction of this modeling procedure provide 
estimates which are somewhat higher than 
conventional modeling^ at higher angles of attack (or 
= 50°, 60*) and slightly lower at lower angles of 
attack (a = 10°, 25°, 30*). For Cj^g E , whose 
results are shown in Figure 12, it was noticed that the 
results of this investigation are consistently lower 


than the previous analysis.^ On the other side the 
results for C N5E in [13] were higher then expected. 
Some reservations about the goodness of the results 
in [13] also existed for and at lower 

angles of attack (a = 10*, 25 ). The results of this 
investigation for C ma and C mfiE , shown in Figures 
13 and 14, reveal more predictable values at these 
angles of attack. 


Conclusions 

Flight data provided by NASA Dryden for 
the longitudinal maneuvers at varying angles of attack 
were used in this preliminary study to determine the 
aerodynamic stability and control derivatives using 
the Maximum Likelihood method for the NASA F/A- 
18 HARV at high angles of attack. The parameter 
identification (PID) code pEst, developed at NASA 
Dryden, was provided by NASA and modified for an 
alternative modeling approach for high angle of attack 
conditions. Preliminary results for key longitudinal 
derivatives using this method were obtained using a 
single set of flight data. A partial proof of the 
suitability of this alternative method for modeling 
high angle of attack conditions is given by more 
meaningful estimates for C ma with respect to classic 
modeling procedures. Additional investigations are 
scheduled using additional flight data. The PID 
process will be extended to the other longitudinal 
derivatives and to the lateral -directional derivatives. 
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Figure 3: Block diagram of the parameter estimation process. 








q (deg/tec) 6 








33 H 





Figure 9: Comparison of time histories for a = 30°. 
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Figure 10: Comparison of time histories for a = 50°. 














Appendix B 


Plots of estimates and time histories 
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